Coherent imaging with pseudo-thermal incoherent light 
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We investigate experimentally fundamental properties of coherent ghost imaging using spatially 
incoherent beams generated from a pseudo-thermal source. A complementarity between the coher- 
ence of the beams and the correlation between them is demonstrated by showing a complementarity 
between ghost diffraction and ordinary diffraction patterns. In order for the ghost imaging scheme 
to work it is therefore crucial to have incoherent beams. The visibility of the information is shown 
for the ghost image to become better as the object size relative to the speckle size is decreased, and 
therefore a remarkable tradeoff between resolution and visibility exists. The experimental conclu- 
sions are backed up by both theory and numerical simulations. 



I. HISTORICAL OVERVIEW AND 
INTRODUCTION 

A decade has passed since the first experimental ob- 
servation of unusual interference fringes in the coinci- 
dence counts of photon pairs [ij, |2(. Signal and idler 
photons produced by parametric down-conversion (PDC) 
were spatially separated and in the signal photon arm a 
double slit was inserted. While no first order interference 
pattern was visible behind the slit, an interference pat- 
tern was observed in the coincidence count by scanning 
the idler photon detector position. To this phenomenon 
it was given the name ghost diffraction. Shortly after 
a ghost image experiment was performed |3(|, showing a 
sharp image of an object placed in the signal arm by reg- 
istering the coincidence counts as a function of the idler 
photon position. 

In the interpretations of the experiments, the quantum 
nature of the source there employed was emphasized, al- 
though the Authors of 3| suggested that "it is possible 
to imagine some type of classical source that could par- 
tially emulate this behavior". Several years passed be- 
fore a systematic theory of ghost imaging (GI) started to 
be developed, and soon a lively debate arose discussing 
the role of entanglement versus classical correlation in 
GI schemes. In the first theoretical papers by the Boston 
group |4( , it was claimed that entanglement was a crucial 
prerequisite for achieving GI, and in particular coherent 
GI: "the distributed quantum-imaging scheme truly re- 
quires entanglement in the source and cannot be achieved 
using a classical source with correlations but without en- 
tanglement. " Soon after, at Rochester University a ghost 
image experiment was performed exploiting the classical 
angular correlation of narrow laser pulses [5| . This fueled 
the debate: which are the features of ghost imaging that 
truly requires entanglement? The debate was continued 
by the paper , where some of us showed that a classi- 
cal GI scheme can indeed produce either the object im- 
age or the object diffraction pattern, but suggested that 
both cannot be produced without making changes to the 
source or the object arm setup. We argued that only 
entangled beams can give both results by only changing 
the setup in the reference arm (the one where the object 
is not present). While by now we know this is not true, 



at that time it was in partial agreement with Q and 
When the Rochester group recently completed the results 
by showing that also the object diffraction pattern can be 
reconstructed using classically correlated beams , they 
had indeed to change the setup (the object location, the 
lens setup as well as the detection protocol). 

Our claim in Q originated from the fact that only en- 
tangled beams can have simultaneously perfect spatial 
correlation in both the near and the far field (in both 
position and momentum of the photons), and no classi- 
cal beams can mimic this l&L In the same spirit, recent 
experimental works El El , brilliantly pointed out a 
momentum-position realization of the Einstein-Podolsky- 
Rosen (EPR) paradox using entangled photon pairs pro- 
duced by PDC. The product of conditional variances in 
momentum and position was there shown to be below the 
EPR bound that limits the correlation of any classical 
(non-entangled) light beam. Based on these results, the 
Authors of pj proposed the same EPR bound as a limit 
for the product of the resolutions of the images formed 
in the near and in the far-field of a given classical source, 
and both EJ argued that in ghost imaging schemes 
entangled photons allow to achieve a better spatial reso- 
lution than any classically correlated beams. 

This was the state of the art, when some of us had 
an idea leading to a ghost imaging protocol with clas- 
sical thermal-like beams. Inspired by the fact that the 
marginal statistics of the signal or idler beam from PDC 
is of thermal nature, we asked ourselves what would have 
been the result of splitting a thermal beam on a macro- 
scopic beam splitter and using the two outgoing beams 
for ghost imaging. Honestly speaking, we expected at 
that time that this would have lead to identify relevant 
differences with quantum entangled beams, where the 
correlation is of microscopic origin. The picture that 
came out was however rather different. The two out- 
put beams of the beam splitter are obviously each a true 
copy (on a classical level) of the input beam: if there is 
a speckle at some position in the input beam, then each 
output beam has also a speckle at the same position. 
Hence the beam splitter has created beams with a strong 
spatial correlation between them, while each beam on its 
own is spatially incoherent. In theoretical works [Til ] we 
showed that this correlation is preserved upon propaga- 
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tion (so it is present both in the near and the far-field 
planes), and that the beams could therefore be used to 
perform GI exactly in the same way as the entangled 
beams from PDC. Actually a very close formal analogy 
was demonstrated between GI with thermal and PDC 
beams, which implied that classically correlated beams 
were able to emulate all the relevant features of quantum 
GI, with the only exception of the visibility [Tlj . 

Thus, we actually had to conclude that what we had 
written in [6j was not wrong but not very correct either. 
What we failed to recognize there is that ghost imaging 
protocols do not need a perfect correlation at all: with 
the imperfect (shot-noise limited) spatial correlation of 
thermal beams both the object image and the object 
diffraction pattern can be reproduced without making 
any changes to the source and only changing the reference 
arm. Moreover the formal analogy between thermal and 
PDC beams suggested that identical performances with 
respect to spatial resolution should be achieved by the 
quantum and classical protocols, provided that the spa- 
tial coherence properties of the two sources were similar. 
This was obviously a controversial result compared to 
what was published until that time 0, 0, and in order 
to be accepted needed an experimental confirmation. We 
recently provided this [12] showing high-resolution ghost 
image and ghost diffraction experiments, performed by 
using a single source of pseudo-thermal speckle light di- 
vided by a beam splitter. As predicted, it was possible to 
pass from the image by only changing the optical setup 
in the reference arm, while leaving the object arm un- 
touched. Moreover, the product of spatial resolutions of 
the ghost image and ghost diffraction experiments was 
there shown to overcome the EPR bound which was pro- 
posed to be achievable only with entangled photons by 
former literature 0, llOj . The origin of the apparent con- 
tradiction with the former literature was there identified, 
by recognizing that the spatial resolution of GI protocols 
do not coincide in general with the conditional variance, 
so that the product of the near and far-field resolution is 
free from any EPR separability bound. 

The idea of using pseudo-thermal light for GI had some 
enthusiastic followers |13l Il4j with proposals for X-ray 
diffraction |15| . some partially converted fans, with ex- 
periments characterizing a pseudo-thermal source of pho- 
ton pairs [T^ and using them for realizing a ghost im- 
age jl7|. and some sceptics 0. The use of pseudo- 
thermal light in GI schemes inspired also a topic which 
became of a certain interest, known as "quantum lithog- 
raphy with classical beams" , or " sub- wavelength inter- 
ference with classical beams" . The quantum version of 
this started with the famous paper by Boto et al. 0] 
claiming that N-photon entangled states could be used 
for improving the resolution of lithography by a factor 
of N. A proof-of-principle experiment using N = 2 in the 
PDC case was provided by [2(j where a halving of the pe- 
riod of the interference fringes was observed in a "ghost 
diffraction" pattern. In Q some of us observed that the 
same effect may be observed when thermal-like beams 



are used, and that in both the entangled and thermal 
case the sub-wavelength interference relies on a simple 
geometrical artifact. We therefore questioned whether 
the Shih experiment really proves Boto's entangled pro- 
tocol. Sub-wavelength interference using thermal beams 
was then theoretically discussed in |2*lj . and experimen- 
tally demonstrated |22| . 

In this paper we continue the experimental investiga- 
tions started in The main result there established 
was that high resolution ghost image and ghost diffrac- 
tion could be achieved with the same classical source, 
with the product of resolutions well behind the EPR 
bound proposed by 0. Here, we shall investigate first 
of all a fundamental complementarity between coherence 
and correlation which exists for ghost imaging schemes. 
Only when the beams are spatially incoherent the corre- 
lation functions allows to retrieve information about the 
object (ghost image or ghost diffraction), while the infor- 
mation disappears as the spatial coherence of the beams 
increases. This is just the opposite of what happens for 
direct detection of the light behind the object, where fully 
coherent information can be obtained only for spatially 
coherent beams. Thus, the spatial incoherence plays an 
essential role for realizing ghost imaging, while instead 
the Hanbury-Brown-Twiss interferometer [23| for deter- 
mining the stellar diameter relies on coherence gained by 
propagation. Secondly, we will investigate visibility and 
signal-to-noise ratio in ghost imaging with thermal light, 
and highlight a tradeoff between visibility and resolution 
when reconstructing the information. 

The paper is organized as follows: In Sec. [H] the ex- 
periment is described, while Sec. II I II introduces the for- 
malism and review the formal analogy between classical 
and quantum ghost imaging. In addition, the relation 
between visibility and signal-to-noise ratio is discussed. 
In Sec. HVI the spatial coherence properties of the beams 
are investigated and experimentally characterized. Sec- 
tion ^focusses on the ghost diffraction setup and shows 
the complementarity between coherence and correlation. 
Section lVll focusses on the ghost image, and discusses vis- 
ibility. In Sec. IVIll numerical results are presented which 
provide a more detailed insight into the results of the ex- 
periment. Finally, the conclusions are drawn in Sec. Willi 



II. DESCRIPTION OF THE EXPERIMENT 

The experimental setup is similar to that of Ref.fl^] 
and is sketched in Fig^ The source of pseudo-thermal 
light is provided by a scattering medium illuminated by 
a laser beam. The medium is a slowly rotating ground 
glass placed in front of a scattering cell containing a 
turbid solution of 3 /im latex spheres. When this is il- 
luminated with a large collimated Nd-Yag laser beam 
(A = 0.532 ^tm, diameter D 10 mm), the stochastic 
interference of the waves emerging from the source pro- 
duces at large distance (z « 600 mm) a time-dependent 
speckle pattern, characterized by a chaotic statistics and 
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by a correlation time r co h on the order of 100 ms (for 
an introduction to laser speckle statistics see e.g.[24|). 
Notice that the ground glass can be used alone to pro- 
duce chaotic speckles, whose correlation time depends on 
the speed of rotation of the ground-glass disk and on the 
laser diameter, as in classical experiments with pseudo- 
thermal light |25l l2q ] . Indeed, in some part of the exper- 
iments described in the following it will be used alone. 
This however presents the problem that the generated 
speckle patterns reproduce themselves after a whole tour 
of the disk, which can be partially avoided by shifting 
laterally the disk at each tour. The turbid solution pro- 
vides an easy way to generate a truly random statistics 
of light, because of the random motion of particles in the 
solution, allowing a huge number of independent patterns 
to be generated and used for statistics. Notice that the 
turbid medium cannot be used alone, because a portion 
of the laser light would not be scattered, thus leaving a 
residual coherent contribution. 
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FIG. 1: Scheme of the setup of the experiment (see text for 
details) . 



At a distance Zq — 400 mm from the thermal source, 
a diaphragm of diameter D — 3 mm selects an angular 
portion of the speckle pattern, allowing the formation 
of an almost collimated speckle beam characterized by 
a huge number (on the order of 10 4 ) of speckles of size 
Air w Xzq/Dq w 21 fim 24]. The speckle beam is sepa- 
rated by the beam splitter (BS) into two "twin" speckle 
beams, that exhibit a high (although classical) level of 
spatial correlation pjj. The two beams emerging from 
the BS have slightly non-collinear propagation directions, 
and illuminate two different non-overlapping portions of 
the charged-coupled-device (CCD) camera. The data are 
acquired with an exposure time (1 — 3 ms) much shorter 
than T co h, allowing the recording of high contrast speckle 
patterns. The frames are taken at a rate of 1-10 Hz, so 
that each data acquisition corresponds to uncorrelated 
speckle patterns. 

In one of the two arms (the object arm 1) an object 
about which we need to extract information is placed. 
The object plane, located at a distance 200 mm from the 
diaphragm, will be taken as the reference plane, and re- 
ferred to as the near-field plane (this is not to be confused 
with the source near-field, as the object plane is in the 
far zone with respect to the source). The optical setup 
of the object arm is kept fixed, and consists of a a single 



lens of focal length F = 80 mm, placed at a distance pi 
after the object and qi = F from the CCD. In this way 
the CCD images the far-field plane with respect to the 
object. 

We shall consider two different setups for the reference 
arm 2. In the ghost- diffraction setup, the reference beam 
passes through the same lens F as the object beam, lo- 
cated at a distance q2 — F from the CCD. In Ref.jT^], 
the spatial cross-correlation function of the reference and 
object arm intensity distributions was calculated, and 
showed a sharp reproduction of the diffraction pattern 
of the object, comparable with the diffraction pattern 
obtained by illuminating the object with the laser light 
(see Sec. EJ. 

In the ghost-image setup, without changing anything 
in the object arm, an additional lens of focal length F' is 
inserted in the reference arm immediately before F. The 
total focal length F2 of the two-lens system is smaller 
than its distance from the CCD, being -r- ~ j, + -jrr. It 
was thus possible to locate the position of the plane con- 
jugate to the detection plane, by temporarily inserting 
the object in the reference arm and determining the po- 
sition that produced a well focussed image on the CCD 
with laser illumination. The object was then translated 
in the object arm. The distances in the reference arm 
approximately obey a thin lens equation of the form 
l/(?>2 — dx) + 1/(72 ~ I/F2 32], providing a magnifica- 
tion factor m — 1.2. In Ref.|12l] the intensity distribu- 
tion of the reference arm was correlated with the total 
photon counts of the object arm showing in this case a 
well-resolved reproduction of the image of the object (see 
also Sec. IVIf) . Thus, the setups allows a high-resolution 
reconstruction of both the image and the diffraction pat- 
tern of the object by using a single source of classical 
light. The passage from the diffraction pattern to the im- 
age is performed by only operating on the optical setup 
of the reference arm, which gives evidence for the the dis- 
tributed character of the correlated imaging with thermal 
light. 



III. FORMAL EQUIVALENCE OF GHOST 
IMAGING WITH THERMAL BEAMS AND THE 
TWO-PHOTON ENTANGLED SOURCE 



The basic theory behind the setup shown in Fig.^has 
been explained in detail in Ref. |ll|. The starting point 
is the input-output relations of the beam splitter 

bi(x) = ta(x) + rv(x) , 62 (^) = ra(x) + tv(x) (1) 

where 61 and &2 are the object and reference beams 
emerging from the BS, t and r are the transmission and 
reflection coefficients of the BS, a is the boson operator of 
the speckle field and v is a vacuum field uncorrelated from 
a. The state of a is a thermal mixture, characterized by 
a Gaussian field statistics, in which any correlation func- 
tion of arbitrary order is expressed via the second order 
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correlation function 



T(x, of) — (a 1 (x)a(x)) 



(2) 



Since the collection time of our measuring apparatus is 
much smaller than the time T co h over which the speckle 
fluctuate, all the beam operators are taken at equal times, 
and the time argument is omitted in the treatment. No- 
tice that we are dealing with classical fields, so that 
the field operator a could be replaced by a stochastic 
c-number field, and the quantum averages by statistical 
averages over independent data acquisitions. However, 
we prefer to keep a quantum formalism in order to out- 
line the parallelism with the quantum entangled beams 
from PDC. 

The fields at the detection planes are given by Cj(xi) = 
J dx'hi(xi, x')bi(x') + Li(xi), where L{ account for pos- 
sible losses in the imaging systems, and h\, ft, 2 are the 
impulse response function describing the optical setups 
in the two arms. The object information is extracted by 
measuring the spatial correlation function of the intensi- 
ties {l\(x\)l2(x 2 )) , where Ii(xi) are operators associated 
to the number of photo counts over the CCD pixel lo- 
cated at Xi in the i-th beam. All the information about 
the object is contained in the correlation function of in- 
tensity fluctuations, which is calculated by subtracting 
the background term (Ji(xi)) (/ 2 (x 2 )) • 



G(xx,x 2 ) = (/ 1 (fi)7 2 (x 2 )) - (Ii(fi))(/ 2 (f 2 )) . 
The main result obtained in was 

G{x u x 2 ) = \rt\ 2 

dx[ / dx' 2 hl(xi,x[)h2(x2,x' 2 )T(x,x') , 



(3) 



(4) 



Equation (0J has to be compared with the analogous re- 
sult obtained for PDC beams [a: 



G P dc(^i,^ 2 ) 
dx[ / dx* 2 h 1 (x 1 ,x* 1 )h2(x 2 ,x 2 )r pdc (x'' 1 ,x' 2 



, (5) 



where 1 and 2 label the signal and idler down-converted 
fields ai, a 2 , and 



r pd c(xi,4) = (ai(^i)a 2 (4)) 



(6) 



is the second order field correlation between the signal 
and idler (also called biphoton amplitude). As already 
outlined in ^]|, ghost imaging with correlated ther- 
mal beams, described by Eq.QJ presents a deep analogy 
(rather than a duality) with ghost imaging with entan- 
gled beams, described by Eq.JSJ): (a) both are coherent 
imaging systems, which is crucial for observing interfer- 
ence from an object, and in particular interference from 
a phase object; (b) both perform similarly if the beams 
have similar spatial coherence properties, that is if F and 
r pc ic have similar properties. They differ in a) the pres- 
ence of h\ at the place of hi, which implies some non 



fundamental geometrical differences in the setups to be 
used and b) the visibility, which can be close to unity 
only in the in the coincidence count regime of PDC. We 
define the visibility of the information as 

, G m ax G max /«\ 

" (hh) max ~ [(h)(h) + G] max - U 



In the thermal case G{x\,x 2 ) 



< 



(/i(fi))</ 2 (x 2 )> so 



ing the signal-to-noise ratio (SNR) associated to a ghost 
imaging scheme (see also [27j)- Intuitively, one expects 



that the visibility is never above i. Conversely, in 
the PDC case it is not difficult to verify that the ra- 
tio Gpdc/ (h) (h) scales as 1 + -X- , where (n) is the mean 
photon number per mode (see e.g. Ref. ll]b). Only in 
the coincidence-count regime, where (n) -C 1, the visibil- 
ity can be close to unity, while bright entangled beams 
with (n) 3> I show a similar visibility as the classical 
beams. However, despite never being above \ in the 
classical case, we have shown that the visibility 

is sufficient to efficiently retrieve information. 

The visibility is an important parameter in determin- 

tiojs: 

that the noise associated to a measurement of I\I 2 is pro- 
portional to (Iil-a), being the statistics of thermal nature. 
This noise also affects the retrieval of the ghost image or 
of the ghost diffraction in a single measurement, because 
this is obtained from by subtracting the background 
term. Hence SNR cx G/(7i/ 2 ), and the visibility defined 
by Eq.J7|) gives an estimate of the signal-to- noise ratio 
of a ghost imaging scheme. This picture is confirmed 
by a more quantitative calculation, not reported here, 
performed by using the standard properties of Gaussian 
statistics. By defining AG = ^(O 2 ) - (O) 2 , with O = 
hh - {h){h), w here G := G(xi,x 2 ), h := h(xi), we 
obtained AG » \/3(/i/ 2 ) 2 + 8G(/i) (I 2 ), where quantum 
corrections have been neglected. If the visibility is small, 
as it is often the case, this reduces to AG « ^(hh), 
and SNRxi ^tt^t- 

Of course, after averaging over N independent mea- 
surements SNR(N) — \Z~N SNR, and if collecting a large 
amount of samples does not represent the problem, any 
ghost image/diffraction can be retrieved after a sufficient 
number of data collections. Hence, if the goal is retriev- 
ing information about a macroscopic, stable object, the 
thermal source represents by far a much better deal than 
the entangled two-photon source. Needless to say, if the 
goal is performing a high sensitivity measurement, or us- 
ing the ghost imaging scheme as a cryptographic scheme 
where information need to be hidden to a third party, 
then the issue of SNR becomes crucial, and the two- 
photon entangled source may turn out to be the only 
proper choice. 



IV. SPATIAL COHERENCE PROPERTIES OF 
THE SPECKLE BEAMS 

Relevant for the ghost image and the ghost diffrac- 
tion schemes are the spatial coherence properties of the 
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speckle beams in the object near-field plane, and in the 
far-field plane with respect to the object. These can be 
investigated by measuring the autocorrelation function 
of the intensities. In any plane it holds a Siegert-likc 
factorization formula for thermal statistics |28, 29] : 

<: mu?) ■■) = (mm*)) + jj im^i 2 , («) 

where M is the degeneracy factor accounting for the num- 
ber of temporal and spatial modes detected. Hence, the 
properties of the field correlation function T can be in- 
ferred from the measurement of the intensity correlation. 
In particular wc will be interested in the field correlation 
function at the object near-field plane T n [x 1 x') 1 and in 
the same function at the far- field plane Tf(£, x'). 

We notice the following equalities, which are a trivial 
consequence of the BS input-output relations Q 



(:h(x)h(x>) : ) = |t| 4 (:/(x)/(f) :) 

= ^(hWhtf)) 

) = |r| 4 (: I(x)I(x') :) 



H$)h(x>) 



(9) 



(10) 



where : : indicates normal ordering and I(x) is the in- 
tensity distribution of the speckle beam in the absence of 
the BS. A part from numerical factors, and from the shot 
noise contribution at x = a?, in a given plane the auto- 
correlation function of each of the two beams coincides 
with the intensity cross correlation of the two beams. 
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FIG. 2: (a)Instantaneous intensity distribution I2 of the 
speckle beam in the near-field plane; (b) Auto-correlation 
function of the intensity in (a). The full line is a Gaussian fit 
of the correlation peak, and the data have been normalized 
to the baseline values. 

Figure shows the instantaneous intensity distribu- 
tion of the reference beam in the setup of Fig^Jwith the 
lens F' inserted, so that the data recorded on the CCD 
are the (demagnified) image of the intensity distribution 
in the near-field plane. A large number of speckles are 
visible with a high contrast, due to the short measure- 
ment time. According to Van-Cittert Zernike theorem 
the size of the speckle here is determined by the inverse 
of the source size (the laser diameter Dq) and and by the 
distance z from the source |24| . Ax n oc Xz/Dq = 32 fim. 



Frame (b) in this figure is the radial autocorrelation func- 
tion (JTDJ, calculated as a function of the distance \x — a?|, 
normalized to the product of the mean intensities. The 
baseline corresponds to the product of the mean intensi- 
ties while the narrow peak located around \x — a?| = 
is proportional to |r„(a?, x')\ 2 , where T n is the second 
order field correlation function at the near-field plane. 
This peak reflects the spatial coherence properties of the 
beams at the object plane. Its width is the near-field co- 
herence length Ax n and gives an estimate of the speckle 
size in this plane Ax n ~ 2ma — 36/im. Notice that 
the peak value is slightly smaller than twice the baseline 
value, giving a degeneracy factor M — 1.7. 

Figure[3]shows the instantaneous intensity distribution 
(a) and the intensity auto-correlation function (b) in the 
far-field plane, measured in the focal plane of the lens F. 
The narrow peak in (b) located around \x— a?| = is now 
proportional to \Tf(x,x')\ 2 , and its width (the far-field 
coherence length) gives an estimate of the speckle size in 
this plane. High-contrast speckles are visible also in the 
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FIG. 3: (a) Instantaneous intensity distribution of the speck- 
les Ii in the far-field plane (b) Auto-correlation function of 
the intensity in (a). The full line is a Gaussian fit of the 
correlation peak, and the data have been normalized to the 
baseline values. 

far-field plane. The Van-Cittert Zernike theorem can be 
again invoked to estimate their expected size, being now 
the source size represented by the diaphragm diameter D, 
Axf oc XF/D w 14 /zm [24(. This is in good agreement 
with the estimation from the the correlation function, 
that gives Axf = 2cr/ = 14.2/im. In this case the peak 
value of the correlation function in frame (b) gives a de- 
generacy factor M = 2.2. This is slightly higher than in 
the near field because Ax{ (the size of the spatial mode) 
is smaller and comparable with the pixel side (6.7/im.) 

Because of the identities in Eqs. (|9"jl. (|10|) . the cross- 
correlation (I1I2) in the near and in the far-field coincides 
with the auto-correlation plotted in Fig[2h and in Fig|3jD, 
respectively. Hence a high degree of mutual spatial cor- 
relation is present in both planes, as a consequence of the 
spatial incoherence of the light produced by our source. 
The more incoherent is the light (the smaller the speck- 
les with respect to the beam size), the more localized is 
the spatial mutual correlation function. The more coher- 
ent is the source (the larger the speckles with respect to 
the beam size), the flatter is the spatial mutual correla- 
tion function. As it will become clear in the next two 
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sections, for highly spatially incoherent light, both the 
ghost diffraction and the ghost image can be retrieved 
with high resolution. Conversely, in the limit of spatially 
coherent light no spatial information about the object can 
be extracted in a ghost imaging scheme, that is from the 
intensity cross-correlation of the two beams as a function 
of the pixel position in the reference beam. 

Summarizing, two aspects of our experiment are cru- 
cial i) the spatial incoherence of light, and ii) a measure- 
ment time <C T co h- Notice that the presence in the near- 
field of a large number of small speckles inside a broad 
beam, implies that the light is incoherent also in the far 
field, because Axf oc l/D, while the far-field diameter of 
the beam oc 1/Ax a . 



V. THE GHOST DIFFRACTION 
EXPERIMENT: COMPLEMENTARITY 
BETWEEN COHERENCE AND CORRELATION 

In this section we focus on the ghost diffraction setup 
(FigH without the lens F'). The object is a double slit, 
consisting of a thin needle of 160 (im diameter inside a 
rectangular aperture 690 fim wide. 

In a first set of measurements the source size is Dq = 
10mm, and the object is illuminated by a large num- 
ber of speckles whose size A.x n — 36^m is much smaller 
than the slit separation. The light is spatially incoher- 
ent as described in the previous section. The results are 
shown in the first row of Fig0] Frame (a) is the instanta- 
neous intensity distribution of the object beam, showing 
a speckled pattern, with no interference fringes from the 
double slit, as expected for incoherent illumination [28| . 
At a closer inspection, the shape of the speckles resembles 



the interference pattern of the double slit, but since these 
speckles move randomly in the transverse plane from shot 
to shot, an average over several shots displays a homo- 
geneous broad spot (Fig^Ib)). Frame (c) is a plot of 
G(xi, X2) as a function of the reference pixel position X2, 
and shows the result of correlating the intensity distri- 
bution in the reference arm with the intensity collected 
from a single fixed pixel in the object arm. Notice that 
at difference to what was done in 12], no spatial aver- 
age [3(j is here employed: this makes the convergence 
rate slower but the scheme is closer to the spirit of ghost 
diffraction in which the information is retrieved by only 
scanning the reference pixel position. The ghost diffrac- 
tion pattern emerges after a few thousands of averages, 
and is well visible after 20000 averages. This is confirmed 
by the data of Fig[SJi which compare the horizontal sec- 
tion of the diffraction pattern from a correlation mea- 
surement to that obtained with laser illumination. In 
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FIG. 5: Horizontal sections of the correlation G(xi,X2) as a 
function of X2, for a fixed x\ (see Fig^t, f). (a) Is the case of 
incoherent light, Do = 10mm; the data are obtained with an 
average over 20000 shots (triangles) and 50000 shots (circles), 
(b) Is the case of partially coherent illumination, Do = 0.1mm 
(20000 shots). The light full line is for comparison the diffrac- 
tion pattern observed with a laser. 
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FIG. 4: Ghost diffraction setup: transition from incoherent 
light to partially coherent light. In the three upper frames 
(a-c) the source size is Do = 10mm, with near-field speckles 
Aa; n = 36/^m. In the three lower frames (d-f) the source size 
is Do = 0.1mm, with Aa; n = 3.2mm. (a)and (d): Instan- 
taneous intensity distribution 7i of the object beam, (b)and 

(e) : Intensity distribution (Ii), averaged over 350 shots (c)and 

(f) : Correlation G(x\,X2) as a function of xi, for a fixed x\, 
averaged over 20000 shots. 
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FIG. 6: Horizontal sections of the average intensity distri- 
bution {Ii(xi)) in the object arm (see Fig^J), e). . (a) Is 
obtained for incoherent light,with Do = 10mm (350 shots), 
while (b) plots the case of partially coherent illumination, 
with Do = 0.1mm (500 shots). The light full line is for com- 
parison the diffraction pattern observed with a laser. 

a second set of measurements the source size is reduced 
to Do = 0.1mm, by inserting a small pinhole after the 
ground glass. As a result, the spatial coherence of the 
light illuminating the object is increased. As the speckle 
size at the diaphragm D is now 3mm , on average 
the object is illuminated by a single speckle of size much 
larger than the slit separation. The results are reported 
in the second row of FigQ] As expected [2^| the interfer- 
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ence fringes are now visible in the instantaneous inten- 
sity distribution of the object beam 1 [frame (d)], and 
become sharper after averaging over some hundreds of 
shots [frame (e)]. Notice that the shape of the interfer- 
ence pattern is now elongated in the vertical direction, 
because the light emerging from the small source is not 
collimated. Horizontal sections of (Ii) , plotted in Fig|5p, 
show a very good agreement with the diffraction pattern 
from laser illumination. Instead, no interference fringes 
at all appear in the correlation function of the intensities 
in the two arms, when plotted as a function of $2 [frame 
(f)]. Notice that in this set of measurements the tur- 
bid medium was removed in order to increase the power. 
This is feasible in this case, because the very small size of 
the source allows a large number of independent patterns 
to be generated in a single tour of the glass disk. 

Figures I4I5I6I evidence a remarkable complementarity 
between the observation of interference fringes in the cor- 
relation function (ghost diffraction), and in the intensity 
distribution of the object beam (ordinary diffraction). It 
also shows the fundamental role played by the spatial in- 
coherence of the source in producing a ghost diffraction 
pattern: the more incoherent is the source, the more the 
two beams are spatially correlated and the more infor- 
mation about the object is available in the ghost diffrac- 
tion pattern. The more coherent is the source, the flat- 
ter is the spatial correlation function of the two beams 
and the less information about the object is contained 
in the ghost diffraction. This is completely analogous to 
the complementarity between the one-photon and two- 
photon interference in Young's double slit experiments 
with photons from a PDC source IP, which was ex- 
plained as a complementarity between coherence and en- 
tanglement. In our case of thermal beams, the comple- 
mentarity is rather between coherence and spatial cor- 
relation, showing that also in this respect the classical 
spatial correlation produced by splitting thermal light 
play the same role as entanglement of PDC photons. 

These results can be easily understood by using the 
formalism developed in Sec lIIII and in particular by in- 
spection of Eq.Qfor the correlation function of the in- 
tensity fluctuations G(xi,X2)- In the limit of spatially 
coherent light the field correlation function T n (xi, X2) be- 
comes constant in space in the region of interest, and the 
two integrals in Eq.(QJ factorize to the product of two 
ordinary imaging schemes, showing the diffraction pat- 
tern of the object only in the object arm 1. As a result, 
by plotting the correlation as a function of X2, no ob- 
ject diffraction pattern can be observed, that is, no ghost 
diffraction occurs. The same observation can be made 
with respect to r p d c (a;i, X2), and G v & c {x\, X2), explain- 
ing thus the analogy between the role of light coherence 
in the PDC and in the thermal case. 

In general, the result of a correlation measure- 
ment is obtained by inserting into Eq. Q the 
propagators that describe the ghost diffraction setup: 
/ii(fi,xi) = {iXF)- 1 e-^ 2l <T(x' 1 ), with T(x) be- 
ing the object transmission function, and h2(x2,x 2 ) = 



{i\F)~ 1 e~ 2 & S2 ' s ' 2 . We get 



G(Xi,X2) oc 



d£T 



2tt 
XF 



. (11) 



where f(q) = J ^e~^ s T(x) is the amplitude of the 
diffraction pattern from the object. The result of a cor- 
relation measurement is a convolution of the diffraction 
pattern amplitude with the second order correlation func- 
tion in the far-field. Hence the far-field coherence length 
determines the spatial resolution in the ghost diffrac- 
tion scheme: the smaller the far-field speckles, the better 
resolved is the pattern. In the limit of speckles much 
smaller than the scale of variation of the diffraction pat- 
tern 



G{xi,x 2 ) 



T 



- v2tt 
X2) XF 



(12) 

Since the amplitude of the object diffraction pattern 
is involved in Ea. HllI) . ghost diffraction of a pure phase 
object can be realized with spatially incoherent pseudo- 
thermal beams, a possibility which was questioned in a 
recent letter [l8|. 



VI. THE GHOST IMAGE: TRADEOFF 
BETWEEN RESOLUTION AND VISIBILITY 

By simply inserting the lens F' in the reference arm 
(see Fig.l), without changing anything in the object arm, 
we now pass to the ghost image. As it was predicted in 
, and experimentally demonstrated in |l2j | , the result 
of cross-correlating the intensities of the reference and 
object arm is now the image of the object, shown e.g. in 
Fig. [7Ji. 

Two issues are important in any imaging scheme: the 
spatial resolution and th e sig nal-to-noise ratio. 

As pointed out in 0, Il2j . the resolution capabilities 
of the ghost image setup are determined by the near field 
coherence length Ax n (the size of the near-field speckles). 
This can be easily understood by inserting the propagator 
h2{x2,x' 2 ) — m5(mx 2 + ^2)1 m t° Eq. Q: 



G(xi, X2) oc 



dfir„(£i, -mx2)T*{x , 1 )e i ^T° 



which shows that the result of a correlation measurement 
in this setup is a convolution of the object transmission 
function with the near-field correlation function r„. In 
the following we shall consider a bucket detection scheme, 
where the reference beam intensity I2 is correlated to 
the total photon counts in the object arm, that is, in 
practice to the sum of photo counts over a proper set of 
pixels. This makes the imaging incoherent [30(, because 
it amounts to measuring 



J dxiG(xi,x 2 ) = J dfj |r n (f 1 , 



-mx2)\ 2 \T(x[)f 



(14) 
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If we take the limit of spatially coherent light, where 
r n (a^,— mx 2 ) can be considered as constant over the 
whole beam size, Both Eqs. l)lH [l .l)14 f show that no infor- 
mation about the object image can be obtained by scan- 
ning X2- Conversely in the limit of spatially incoherent 
light where the speckle size is much smaller than the scale 
of variation of the object, T n (x' 1 , —771x2) ~ <5(x*i + 
and both Eqs. |T3 )l .lfH jl converge to \T (— 77iX2)| 2 xconst.. 

Concerning the signal-to-noise ratio, the discussion in 
Sec. IIIII pointed out that it is determined by the image 
visibility. We have studied the visibility of the ghost im- 
age of a double slit in a sequence of measurements where 
the vertical size of the apertures was progressively re- 
duced, while leaving unchanged their horizontal size and 
separation. This is shown in Fig|JJi> where all the frames 
display the correlation function measured in a bucket de- 
tection scheme as a function of the reference pixel posi- 
tion x 2 . Despite the fact that all the frames have been 
obtained with the same number of averages N = 10000, 
the sequence display a remarkable enhancement of the 
visibility as the object area is reduced. This enhancement 





FIG. 7: Reconstruction of the object image via correlation 
measurements (FigQ with the lens F 1 inserted), (a) Cross- 
correlation of the intensity distribution of the reference arm 
with the total photo counts in the object arm, as a function 
of X2- (statistics over 10000 CCD frames). In the sequence 
of frames the object area is progressively reduced, and cor- 
respondingly a enhancement of the visibility can be observed 
(b) Horizontal sections of the images in (a), with the correla- 
tion normalized to (I1I2), so that the vertical scales gives the 
visibility. 

is more clearly visible in the horizontal sections plotted in 
Fig. |7p, where in each point the correlation function has 
been normalized to (I1I2), so that the numbers on the 
vertical axis give directly the visibility [see Eq. Q]. We 
notice that the visibility increases as the object area de- 
creases, and correspondingly the the signal-to-noise-ratio 
increases, as expected. 

In order to understand the origin of the behaviour 
shown in Fig0 we need first to consider Eq.JTlJ, that 
gives the correlation function in the bucket detection 
scheme. Let us assume that the object simply trans- 
mits the light over a region of area S bj and stops it 
elsewhere. By assuming that the coherence length Ax n 
is smaller than the object features, as it is necessary for 
the object to be resolved, the integrand at r.h.s of ifH|) 



can be non-zero for x' x in a region located around x*2, of 
area A co h, where A co h is the coherence area oc Ax^. Thus 
the correlation scales as the coherence area. 



dxiG(xi,x 2 ) oc A C oh 



(15) 



Conversely, it is not difficult to see that in the bucket 
detection scheme 



dx^hixi)) = / dx[ |T(fi)r (J n (£i)> oc A obj , (16) 



where (^(x^)) is the average intensity distribution of the 
light illuminating the object, that can be taken as roughly 
uniform on the object area (the speckles average to a 
broad uniform light spot, as shown in Fig. 0J>). Hence 
the ratio of the correlation to the background scales as: 



/ dxiG(x*i,x 2 ) 



oc 



Jdxi(Ii(xi))(7 2 (x2)) A obj 



(17) 



This formula is reminiscent of role of the mode degen- 
eracy in Eq.©, and indeed it reflects the fact that in 
a bucket detection scheme the number of spatial modes 
detected is proportional to A bj/A co h, which represents 
a degeneracy factor that reduces the visibility of the cor- 
relation with respect to the background. The ratio in 
Eq. H17|) is usually small, so that the visibility of the 
ghost image roughly coincides with it. Thus the visi- 
bility roughly scales as the ratio of the coherence area 
to the object transmissive area: the larger are the object 
dimensions with respect to the speckles, that is, the more 
incoherent is the light illuminating the object, the worse 
is the visibility of the ghost image. This is confirmed by 
the plot in Fig. [5J showing how the visibility increases 
with the inverse of the object area. This rather counter- 
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FIG. 8: Visibility of the ghost image as a function of the 
inverse of the transmissive area of the object, showing an 
increase of the visibility by reducing the object area. 

intuitive result also implies that a better resolution can 
be achieved only at the expenses of the visibility, since 
the resolution is determined by the speckle size. Hence, 
complex images that need small speckles to be resolved in 
their details have a lower signal-to-noise ratio than simple 
images which can be resolved with relatively large speck- 
les (see also ^3 f° r a similar conclusion). This, however, 
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does not prevent from retrieving more complex images 
(see e.g. Fig. |5J), provided that a larger number of data 
acquisitions are performed. 






FIG. 9: Ghost image of the number 4, retrieved from the 
correlation function after averaging over 30000 shots. 



FIG. 10: Comparing a two dimensional numerical simulation 
of the experiment, by showing the correlation of intensity fluc- 
tuations normalized to the product of the beam intensities, 
(a) The ghost diffraction case G(xi, X2)/{h {xi)) (hix'i)) ■ (b) 
The ghost image case J Ax\G(x\,x2) / {I2{x2)) J dxi(Ix(xi)). 
In both cases the numerics and experiment are real units, 
and are as reference compared to the data obtained by coher- 
ent laser illumination of the object. The averages are done 
over 2 ■ 10 4 realizations. In the numerics Ax n = 34 fim and 
Axf = 12 fim. 



VII. NUMERICAL RESULTS 



In this section we use a numerical model for simulating 
the speckle pattern created by the ground glass and the 
turbid medium, as to support the results of the previ- 
ous section. The thermal field is created by generating a 
noisy field with huge phase fluctuations. We then mul- 
tiply the noisy field with a Gaussian profile and a sub- 
sequent Fourier transformation gives what corresponds 
to the near field; the width of the Gaussian then con- 
trols the near-field speckle size. The far-field speckle size 
is controlled independently as in the experiment: after 
generating the near field, a diaphragm of diameter D 
is introduced, beyond which only vacuum fluctuations 
appear; D then controls the far-field speckle size. The 
speckle field transmitted by the diaphragm is then im- 
pinged on a 50/50 BS, with vacuum fluctuations entering 
the unused port, giving the two desired correlated beams. 
We neglect the temporal statistics, since we assume that 
the short measurement time of the experiment provides a 
speckle pattern frozen in time. We should finally mention 
that a Wigner formalism is used to describe the quantized 
fields, as described in detail in 

Initially, let us briefly show that the numerical simula- 
tions are able to describe very precisely the experiment. 
In Fig. ^3 is shown the results of two-dimensional (2D) 
simulations with all parameters kept as close as possible 
to the experiment. This includes near-field and far-field 
speckle sizes, object and aperture sizes, as well as num- 
ber of realizations. Both the ghost diffraction pattern 
[Fig. ITUT a)] and the ghost image [Fig. [TOT b)] show very 
good agreement with the experimental recorded data (us- 
ing the small-speckle setup of Sees. lIVl IVI|) . We stress 
that this comparison is not in arbitrary units. 

In Sec. E] we showed experimentally the behaviour 
of the system when using either coherent or incoher- 
ent beams to investigate the diffraction properties of an 
object. In order to investigate better the actual transi- 
tion from coherent to incoherent illumination of the ob- 



ject, we have carried out numerical simulations that are 
presented in Fig. ^2 We have there kept D = 3 mm 
and then for each simulation changed Ax n . The simula- 
tion only includes one spatial direction (ID), and there- 
fore the coherence properties of the beam is governed by 
the ratio Aa; n /L bj, where L bj is the ID equivalent of 
A bj. Thus, the smaller Acc n /Z/ bj the more incoherent 
is the beam impinging on the object. For small speck- 
les (Ax n /L bj <C 1) the beams are spatially incoher- 
ent, implying a strong spatial correlations between the 
beams: the ghost diffraction is observed in the correlation 
[Fig. lllf a)]. In contrast, no diffraction pattern can be ob- 
served directly in the object beam [Fig. Illf b)]. As Ax n 
is increased by generating bigger speckles the beams be- 
come more and more spatially coherent (Ax n /L bj — 1): 
the ghost diffraction disappears gradually [Fig. Illf a)], 
while the diffraction pattern starts to appear from the 
direct observation of the object beam [Fig. Illf b)]. Fig- 
ure ^^c) shows what happens for the ghost image dur- 
ing this transition: the incoherence for small Ax n implies 
that a ghost image of the object can be observed, and it 
disappears gradually while increasing the coherence. 

We saw in Sec. lVll that the visibility of the ghost image 
became better as the object area was reduced, cf. Fig. |S] 
To investigate this phenomenon in general we show in 
Fig.ll2lhow the object size|33j effects the visibility of the 
information. The trend we saw in the experiment is re- 
peated in the numerics: in Fig. I12f a') the ghost image 
visibility increases as the object size decreases because 
fewer modes are transmitted. In Fig. I12f b~) the simu- 
lation is repeated for the ID case with a similar result. 
However, since in the ID case much fewer modes are 
transmitted by the object the visibility is much higher. 
We have also in Fig.^]plotted the visibility of the ghost 
diffraction fringes, and we observe that - in contrast to 
the ghost image case - the visibility decreases as the 
object size is decreased. This is is to be expected be- 
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FIG. 11: ID numerical simulation of the experiment show- 
ing the transition from incoherent to coherent illumination 
of the object, realized by changing the near-field speckle 
size Ax n - (a) shows the normalized correlation of intensity 
fluctuations in the ghost diffraction case, while (b) shows 
the normalized (Ii(xi)} as observed directly in the object 
arm. (c) shows the correlation of intensity fluctuations 
in the ghost image case, normalized to the beam intensi- 
ties J dxiG(xi, X2)/(h{x2)} J dxi(Ji(xi)}. The averages are 
done over 10 5 realizations. The object mimics the experimen- 
tal one, implying L ob j = 530 fim. Axi = 12 /im. 



cause for the diffraction pattern the modes transmitted 
by the object interfere coherently so G(x\,X2) oc 
(for the 2D case). In contrast, for the mean intensity the 
modes interfere incoherently so (Ii(xi)} oc A y,y Thus 
G(xi, x 2 ) / (h(xi)) (I 2 {x 2 )} oc A obi : the bigger the object 
the better the visibility of the information. We also note 
that there is basically no difference between the ID and 
2D results for the ghost diffraction visibility. This is be- 
cause a similar argument can be done for the ID case 
showing G(a;i,a;2)/(/i(a;i)>(/2(a;2)) oc L ohy Finally, we 
have checked that if the far-field speckle size Axf is varied 
and all other parameters are kept fixed, then the visibil- 
ity of the diffraction fringes increases as Axf is increased: 
again a tradeoff between resolution and visibility is found. 




FIG. 12: Numerical simulations of the experiment showing 
how the object size effects the visibility V . We kept the 
speckle sizes constant but varied the width of the two slits, 
(a) is the 2D case showing V as function of the speckle area 
relative to the object area. Note that the ghost image visibil- 
ity has been multiplied by a factor of 30, and that the object 
length perpendicular to the slits was kept constant, (b) is the 
ID case, showing V as function of the speckle size relative to 
the object length. Ax n = 34 /im and Axf = 12 /im. 



VIII. CONCLUSIONS 

The experimental results reported in this paper con- 
firms that correlated imaging can be performed with a 
classical thermal source. A remarkable complementarity 
between spatial coherence and correlation is predicted 
and confirmed by experiments and numerical simula- 
tions. By changing the coherence of the speckle field at 
the object plane from incoherent to coherent (measured 
relative to the object dimensions), the object diffraction 
pattern reconstructed from correlations disappears while 
it appears in the far field intensity distribution measured 
in the object arm. We also discussed from a quantita- 
tive point of view the issue of the visibility of the cor- 
related imaging scheme. We showed that the visibility 
of the object image is proportional to the ratio between 
the object area and the field coherence area at the ob- 
ject plane. This means that a tradeoff between resolution 
and visibility exists: a better visibility can be obtained 
only at the expense of a lower resolution and vice versa. 
However, the experiment clearly shows that a fairly good 
resolution can be achieved since the problem of low visi- 
bility can be overcome by performing a sufficiently large 
number of averages. 
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